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Context. Because of its large angular size and proximity to the Milky Way, NGC 253, an archetypal starburst galaxy, provides an 
excellent laboratory to study the intricacies of this intense episode of star formation. 

Aims. We aim to characterize the excitation mechanisms driving the emission in NGC 253. Specifically we aim to distinguish between 
shock excitation and UV excitation as the dominant driving mechanism, using Bry, H2 and [Fell] as diagnostic emission line tracers. 
Methods. Using SINFONI observations, we create linemaps of Bry, [Fell]i ^4, and all detected H2 transitions. By using symmetry 
arguments of the gas and stellar gas velocity field, we find a kinematic center in agreement with previous determinations. The ratio of 
the 2-1 S(l) to 1-0 S(l) H2 transitions can be used as a diagnostic to discriminate between shock and fluorescent excitation. 
Results. Using the 1-0 S(l)/2-l S(l) line ratio as well as several other H2 line ratios and the morphological comparison between 
H2 and Bry and [Fell], we find that excitation from UV photons is the dominant excitation mechanisms throughout NGC 253. We 
employ a diagnostic energy level diagram to quantitatively differentiate between mechanisms. We compare the observed energy level 
diagrams to PDR and shock models and find that in most regions and over the galaxy as a whole, fluorescent excitation is the dominant 
mechanism exciting the H2 gas. We also place an upper limit of the percentage of shock excited H2 at 29%. 

Conclusions. We find that UV radiation is the dominant excitation mechanism for the H2 emission. The H2 emission does not correlate 
well with Bry but closely traces the PAH emission, showing that not only is H2 fluorescently excited, but it is predominately excited 
by slightly lower mass stars than O stars which excite Bry, such as B stars. 

Key words. Keywords should be given 



1. Introduction 

Understanding the excitation mechanisms and drivers in nearby 
starburst galaxies is critical to understanding their role in 
galaxy evolution. One of the best laboratories to st udy this 
is NG C 253. This galaxy is a nearby (3.5 Mpc Reko la et alJ 
(I2005h 1. edge-on infrared bright galaxy that is part of the 
Sculptor group. It is classified as a late-type barred spiral 
galaxy (SAB(s)c) with star formation confined to the inner 0.5 
kpc. The star formation rate of the starburst is 1.4-9.5 
yr"' dGrimm e t al. 2003). This intense episode of star forma- 
tion i s driven by a 7 kpc b ar, which funnels gas into the nu- 
cleus (lEng elbracht et a lJl998i) . There is a starburst driven super- 
wind extending perpen dicular to the disk dP emoulin & Burbidge 
ll970';'Ulrich 1978; Fabbiano & Trinchierin984tlMcCarthy et alJ 
11987: WestmoQuette et al. 201 1). 

The inner region of star for mation is highly obscured 
with an Ay of ~ 5 - 18 mag (" Kornei & McCradyl l2009t 
[pernand ez-Ontiyeros et alj r2009). Many dense, luminous, dust 
enshro uded young stellar clusters have been found in the nu- 
cleus jFo^eset^!lll99lj20^ In addition to the stellar clus- 



ters, AIonso-Herrero et alJ (l2003h and lUlvestad & Antonuccil 



(Il997ir have observed many compact radio sources in [Fell] 
and radio resp ectively, representing young supernova remnants. 
iTurner & Hoi d 19851) observed over 60 compact radio sources, 
the brightest of which is named TH2. A high resolution (1"- 
2") study of the molecular gas in NGC 253 was performed by 
ISaka moto et al.l (1201 ih . They found 5 warm, dense clumps of 
molecular gas, which are all coincident with radio sources. They 
suggest that these clumps are natal molecular cloud complexes 
harboring massive star formation. The K band continuum peak. 



which was originally defined as the nucleus by iForbes et al.l 
(1199 Ih . is the location of a super star cluster. This s uper star clus- 
ter wa s observed with Hubble Space Telescope byl Watson et al.l 
(119961) and is consistent with a single stellar population. In the 
mid infrared, two bright peaks have been observed in the cen- 
tral regions. The brightest of the two is coincident with the K 
band contin uum peak and super star clu ster and has no radio 
counterpart dUlvest ad & A ntonuccil 1 19971) . The second, 3.0" to 
the northe ast is coin cident to the radio peak TH2. A more recent 
study by Ulvestad & Antonuccil d 19971) suggests that TH2 is the 
true nucleus of NGC 253 and represents either a low-luminosity 
AGN or a compact supernova remnant. 

Although T H2 is generally considered to be the nucleus, a 
recent study by iMiiller-Sanchez et al.l d201Q) provides new es- 
timates for the kinematic center based on the kinematics of 
the stellar component. Using SINFONI on the VLT, they were 
able to achieve very high angular resolution and derive the stel- 
lar velocity field. They found the kinematic center to be offset 
from TH2 only by Ax=0".6 and Ay=0".4 (i.e. by 12 pc) in the 
southwest direction. They propose an altern ative position for th e 
kinematic center as the X-ray source X-1 dWeaver et al.ll2002h . 
Although NGC 253 has been the target for many detailed kine- 
matic studies, there is still much uncertainty about the precise 
location of the true kinematic center. 

The subject of e xcitation mechanism s has been extensively 
studied in NGC 253. lMartfn et all d2006l) finds that the chemistry 
and heating of NGC 253 is dominated by large scale, low veloc- 
ity shocks. Presence of shocked molecular material is indeed ev- 
ident through the prese nce of widespread SiO emission through- 
out the nuclear regions dGarcfa-Burillo et ani2000l) . lMartm et al.l 
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(l2009bl) suggest that, although NGC 253 is dominated by shock 
chemistry, PDRs play a crucial role in the chemistry, since there 
are very high abundances of PDR tracing molecules, namely 
HCO+, CO+. In addition, they find that although NGC 253 was 
thought to be at an earlier stage of evolution than M 82, a proto- 
typical starburst with strong PDR characteristics, the molecular 
clouds are larger and have a higher column density, pointing to 
a later stage of evolution, such as seen in M 82. Martin et al, 
( l2009ah also find that the UV heating in NGC 253 is more simi- 
lar to M 82 than o ther starbursts, yet stil l not the dominant exci- 
tation mechanism. i Harrison et alj(ll998h used the ortho- to para- 
ratio of H2, to argue that m ost of the H2 is excited by PDRs 
instead of shocks. In addition. Ilsrael et aP (ll989h : lLuhman et alJ 
11991 show in NGC 604 and Orion A respectively that although 
in small -beam apertures it may seem as if high surface brightness 
H2 emission from shocks is dominant, when the beam size in- 
creases, it is clear that fluorescent excitation dominates on large 
scales. 

In this paper, we investigate the excitation mechanisms for 
the molecular gas in the nuclear region of NGC 253. We use the 
following methods to differentiate between shock dominated and 
PDR dominated excitation on a 2 pc scale: 

- Calculate the 2-1 S(l)/l-0 (SI) line ratio at each pixel posi- 
tion 

- Compare morphology of Bry, [Fell], and PAHs to that of the 
H2 

- Relate observed H2 excitation diagrams to those generated 
by shock and PDR models. 

In Section |2] we will discuss the observations, while in 
Section [3] we present the derived spectra and linemaps. In 
Section |4] we will investigate and define the kinematic center. 
Then, in Section |5] we investigate the nature of the dominant 
excitation mechanisms throughout the galaxy, focusing on the 
dense clumps. In Section|6] we compare our excitation diagrams 
to shock and PDR models spanning the parameter space. We 
summarize our conclusions and their implications in Section |7] 

2. Observations 

2.1. SINFONI Observations 

All observations were made with the Spectrograph for INtegral 
Field Observations in the Near-Infrared (SINFONI) at the ESO 
VLT. SINFONI provides spatial and spectral data in the form of 
data cubes in J, H, and K bands. The SINFONI instrument is 
mounted at the Cassegrain focus of the Unit Telescope 4 at the 
Very Large Telescope (VLT). 

We observed in the H, and K bands using a spatial pixel scale 
of 0.25" corresponding to a field of view of 8"x8" per frame and 
a spectral resolution of 2000, 3000 and 4000 respectively, which 
corresponds to a velocity resolution of 149.8, 99.9 and 74.9 
km/s. All science observations were taken in the ABA' nodding 
mode (300s of object, 300s of sky, 300s of object), where A' 
is slightly off'set from A. The object exposures are aligned and 
averaged during the reconstruction of the data cube. 

The observations of NGC 253 were made in visitor mode on 
August 28th, 2005. In order to capture the full extent of the H2 
emission, consecutive frames were taken in the K band moving 
further away from the center, along the disk until H2 was no 
longer detected. This resulted in 6 separate pointings. Since there 
are also H2 transitions in the H band, a similar strategy was used, 
resulting in 4 separate pointings. 



We used the standard reduction techniques of the SINFONI 
pipeline on all observations, including corrections for flat field, 
dark current, nonlinearity of pixels, distortion, and wavelength 
calibration. We obtained the flux calibration and atmospheric 
corrections from observations of a standard star, namely HR 
2058 in the H band and HD 20001 in th e K band. This is the 
same dataset' Miiller-Sanchez et aT ] (120101) used to determine the 
kinematic center. 

2.2. ISAAC Observations 

We use the ISAAC observations of NGC253 in the 3.21 and 3.28 
jum filters described by Tacconi-Garman et al. (2005). These ob- 
servations were uncalibrated. We established a flux calibration 
for these two images using the ISO-SWS spectra by Sturm et al 
(2000), taking into account the filter profiles of the two ISAAC 
filters and the precise location and orientations of the ISO-SWS 
aperture. Inspection of the spectrum by Sturm et al (2000) shows 
that the 3.21 yum filter contains only continuum emission, while 
the 3.28 /im filter contains continuum+PAH emission. The peak 
of the PAH emission is at 3.3 //m, which is fully included in the 
3.28 jjm filter The spectral slope is small enough that there is 
negligible change in continuum level between 3.21 and 3.28 //m 
and therefore a PAH emission image was created by subtracting 
the calibrated 3.21 /im image from the calibrated 3.28 yum image. 

3. Results 

The SINFONI datacubes allow us to construct continuum maps, 
line maps and velocity maps of emission lines detected in these 
bands. The K band continuum map and H2 1-0 S(l) emission 
line map are shown in Figure [T] and the Bry and [Fell] 1 54 emis- 
sion line maps are shown in Figure |2] All line maps are cen- 
tered on the kinematic center proposed by Miiller-Sanchez et aTj 
(l2010h . We have chosen 8 regions defined by the most domi- 
nant H2 peaks and these regions are overplotted by white rectan- 
gles in each map. There is an additional region centered on the 
kinematic center All regions are 12x12 pixels or 1.5x1.5 arc- 
seconds. These regions were selected based on H2 flux. There 
is a relatively bright H2 emission clump to the north west of 
Region 1, which was not selected for analysis since it shows 
similar emission features to Region 1 . 

The K band continuum map reveals the inclined disk of the 
galaxy, tracing the older stellar population. The disk appears 
irregular with clear clumps of high luminosity. The kinematic 
center (Region 3) shows a very high level of continuum emis- 
sion. In addition, superimposed on the bright continuum emis- 
sion around the kinematic center, there is a bright flux peak 
at the location of th e super star cluster (Region 2) studied by 
IWatson et al.l(ll996h . 

On the other hand, the H2 1-0 S(l) map shows a completely 
different morphology. The H2 extends along the entire disk with 
its brightest clump located several arcseconds east of the kine- 
matic center (Region 3). However, the H2 is even more inhomo- 
geneous than the K band continuum. There are clear H2 bright 
clumps ranging out into the furthest regions. The white rectan- 
gles highlight the brightest emission clumps. There is also signif- 
icant H2 emission originating from the super star cluster (Region 
2), but not from the center. 

Similarly, comparison of the H2 and the Bry map reveals 
little correlation. The Bry flux peak is coincident with the su- 
per star cluster, but shows no clear flux increase near Region 4, 
where H2 is brightest, nor at the center. The Bry is also gener- 
ally less extended than the H2, with most of the emission in the 
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Fig. 1. K band continuum map and H2 1-0 S(l) line map of NGC 253. The offset from the kinematic center at position (a, 5)2000 = 
(00H7"'33M4, -25°17'17".52) is given on each axis. The contours represent 13 equally spaced levels with a maximum flux of 
1.3x10"" erg s^'cm^^/^m^'arcsec"^ and5.4x 10"'^ erg s"'cm"-arcsec"^,for the Kband continuum and H2 respectively. The white 
rectangular regions are H2 flux peaks, with the exception of Region 3, which is centered on the kinematic center. The blue '+3 
represents the kinematic center, the blue triangle represents X-1 and the blue diamond represents TH2. 
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Fig. 2. Bry and [Fell] 1 64 emission line map of NGC 253. The offset from the kinematic center at position (a, 5)2000 = 
(00H7"'33M4, -25°17'17".52) is given on each axis. The contours represent 13 equally spaced levels with a maximum flux of 
4.5x10"'"* and 2.8x10"''* erg s"'cm"^arcsec"^ for Bry and [Fell] respectively. The white rectangular regions are areas of concen- 
tpated H2 flux, with the exception of Region 3, which is centered on the kinematic center. 
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southwestern portion of the galaxy, centered on the super star 
cluster. 

The [Fell] 1.64 map also reveals emission peaking near, but 
not coincident with. Region 4. In the northeastern regions, the 
[Fell] emission is also quite clumpy, with some of the clumps 
coinciding with the H2, such as Regions 5, 7 and 8. Comparing 
the [Fell] to the Bry emission, we see in both cases, as well as 
in the K continuum, an elongated arm-like feature coming from 
the center, curving out towards the northeast corner along the 
disk, bordering on the northern edge of Region 5. However the 
[Fell] feature is spatially offset toward s the south wi t h resp ect to 
the Bry feature, as noted originally bv lForbes et alj (11993 ). The 
H2 and continuum emission in this region traces neither Bry or 
[Fell]. 

Figure [3] displays the rest-frame K band spectra, integrated 
over the regions shown on the H2 map (Figure [T]i. The main 
emission lines are highlighted. Specifically, we see many H2 
transitions as well as Bry and Hel. The highest energy H2 transi- 
tion detected is the 3-2 S(3) transition. This is located at 2.2008 
yum and only marginally visible in Region 8 of Figure[3] although 
we have 2cr detections of it in Regions 1, 2, and 3. This is an im- 
portant detection since this transition can only be excited by high 
energy photons, and therefore points towards fluorescent excita- 
tion. In Table [T] we present the integrated fluxes and fluxes in 
each region of all detected H2 lines with error bars. These error 
bars were derived by calculating the flux for a range of different 
continuum baselines and finding the standard deviation of these 
values. We also attempted to extract H2 line fluxes from transi- 
tions in the H band, yet the stellar absorption features overpow- 
ered the emission lines and we were unable to get reliable flux 
measurements. 

Using the SINFONI data cubes, we can extract both a gas 
and stellar velocity map. Figure |4] displays the stellar veloc- 
ity field with the H2 molecular gas velocity field in contours. 
The stellar velocity field was calculated by fitting the three stel- 
lar CO photospheric absorption features. Fitting with a set of 6 
template stars (K3V, MOIII, MOV, M4y , M5III and M5II), we 
used the Cap pellari & E msellem' (2004) Penalised Pixel Fitting 
(PPXF) package. All the template stars were observed with 
SINFONI on the same settings as the NGC 253 observations. 
In a ddition, t he pixels wer e binned using a Voronoi tessellation 
(ICap pellari & Copinll2003h . with a signal-to-noise of 50 per bin. 
These binned spectra were then fit with the best combination of 
the 6 template stellar spectra varying velocity and dispersion. 
Then, the systemic velocity of the galaxy, 243 km/s, was sub- 
tracted from the velocity field. 



4. The Kinematic Center of NGC 253 

The stellar velocity field and the H2 velocity field show many 
general similarities, especially in the central region. Specifically, 
both the stellar and gas velocity fields show an 'S' -shaped mor- 
phology in the velocity contours in the km/s range. This mor- 
phology is typical of galaxies with a bar potential, and supports 
the claim of a bar potential in NGC 253. 

Since the H2 velocity field is much cleaner than the stel- 
lar velocity field and they share a similar morphology, we will 
use the H2 velocity to estimate the kinematic center of NGC 
253. In order to estimate this position, we create a mirrored 
reflection around a set of central points that could be the cen- 
ter By subtracting this mirrored map from the original map, 
we can find which pixel provides the most symmetric profile. 
Using this method, we can determine the center pixel, or the 



pixel with the most symmetry. To determine the absolute coor- 
dinates of this pixel, we employ the 2MASS K band image of 
NGC 253. We anchor the SINFONI observations to the 2MASS 
observations by finding the position of the K-band peak. We then 
use the ofi'set of the kinematic center to the K band center to 
determine the a bsolute coordinates, (a g)2Qo o=(00''47'"33 '.084,- 
25°17'18".42'). lMuller-Sanchez et al.l (12010 ) find the kinematic 
center at (o-, (5)2ooo=(OOM7'"33M4,-25°I7' 17".52) with a mean 
3cr error of r=1.2". Thus, our center is within this error and in 
good agreement with their result. The diff'erence between the two 
results is due to the method. Muller-Sanchez used an inclined 
disk model to find their kinematic center, however the symmetry 
argument focuses on finding the point of highest symmetry. This 
method is compromised due to two symmetric circular regions 
in the velocity field, the top circle is located within the red error 
circle and the symmetric counterpart is just below it and appears 
green. Our method finds the center to be just between these two 
circles, at the radius of the Muller-Sanchez error radius, repre- 
sented as a square. Due to this symmetric point, we consider the 
Muller-Sanchez center to be the true kinematic center, which is 
also clear from the inflection point of the bar profile 'S' shape 
seen in the H2 velocity field. Since the H2 inflection point co- 
incides perfectly with the Muller-Sanchez center, it proves that 
there is no true ofifest between the H2 and CO measured cen- 
ters. We therefore independently confirm the kinematic center 
proposed by Muller-Sanchez et al. (2010). In FigureEl the H2 I- 
S(l) velocity map is presented with the lMiiller-Sanchez et al.l 
( 20 1 0) 3 cr center circled in red and their proposed kinematic cen- 
ter overlaid with a white plus sign. For completeness, the radio 
peak TH2 and the x-ray peak X-1 are both overplotted as a dia- 
mond and triangle respectively. 

5. Gas Excitation l\/lechanisms in NGC 253 

In general, the morphologies of the K band continuum, H2, Bry, 
and [Fell] 1 64 emission are dissimilar. None of the peaks are co- 
incident with each other with the exception of the K band and 
Bry peaks at the super star cluster There are certain regions 
where the H2 matches the Bry, while some H2 clumps match 
the [Fell] 1 64 emission. This suggests that the H2, Bry and [Fell] 
are being excited through different physical processes. 

The Bry peak, tracing massive, young star formation, is co- 
incident with the K band continuum peak, at the location of the 
super star cluster. In general, the K band continuum traces the 
older medium and low mass populations of stars in contrast to 
the youngest, most massive stars powering the Bry emission. 
There is extended emission in the northeastern direction of the 
Bry peak, but no other strong emission regions. 

[Fell] is another strong NIR diagnostic line, which traces 
shocks. In interstellar space, iron atoms are typically fully locked 
into dust grains. Shock fronts associated with supernova rem- 
nants (SNRs) may cause very efficient grain destruction through 
thermal sputtering dJones et alJll9 94). This releases the iron into 
the gas-p hase where it is sin gly ionized by the interstellar radia- 
tion field (iMouri et alJl2000l) . In the extended post shock region, 
Fe^ is excited by electron collisions, rendering it a strong diag- 
nostic shock tracer. The [Fell] peak is located slightly northeast 
of the H2 peak (Region 4). It shows diffuse emission consistently 
along the disk, but does not show a ny peak at Region 2, the lo- 
cation of the super star cluster. From [Rosenberg et al.l (|20I2|) we 
can calculate the supernova rate (SNrate) from the Fell lumi- 
nosity. In Table |2] the calculated SNrates for each region and 
the whole galaxy is given. The small amount of [Fell] flux and 
correspondingly low supernova rate suggests that the super star 
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f'ig. 3. K band spectra of H2 bright clumps. The integrated K band spectra, in the rest frame, integrated over the regions shown in 
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Table 1. Integrated fluxes for selected regions and over the full map of all detected H2 transitions. For the regions the units of flux 
are 10 erg s 'cm^^. Only detections above 3cr are reported with the exception of 1-0 Q(4), where we report all fluxes over 2cr. 
The errors on the integrated fluxes represent 1 standard deviation. Each region is 12x12 pixels, or 1.5x1.5 arcseconds. 
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Fig. 4. Stellar velocity field (left) with H2 1-0 S(l) velocity contours and H2 1-0 S(l) velocity map with contours (right). The contour 
levels are -75 ,-50 ,-30 ,-10 ,0 ,10 ,30 ,50 ,75 km/s. The color bar describes the H2 a nd stellar velocity in km/s. On the right panel, 
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circle, and the kinematic center. The diamond is TH2, the radio peak. The triangle is X-1, the X-ray peak. The square is the point of 
highest symmetry. 



cluster in Region 2 is a young cluster. 'Fern andez-Ontiveros et alJ 
(12009 ) suggests that Region 2 is a starburst of 6 Myrs, which is 
also consistent with our age calculations. Our integrated super- 
nova rate in the inner 300 pc of the galaxy, is 0.2 yr which is 
wefl matched to previous supernova rate measurement of 0. 1 -0.3 
(iRieke et al.l[T98a [19881: luivestad & Antonuccill 19971) . 

Unlike Bry and [Fell], the mechanism which excites the H2 
gas is not as clear. The NIR is rich with ro-vibrational H2 lines, 
which can either be excited thermally, through shocks, or by flu- 
orescence, through UV photons from O and B stars. In the case 
of shock excited H2 emission, the gas is thermalized and the en- 
ergy levels are populated in a "bottom-up" manner. The result- 
ing gas temperature is around 1000 K and thus there is little H2 
emission from the v > 3 states, which have temperatures above 
15,000 K. However, in the case of excitation by UV photons, 
the H2 molecule absorbs a highly energetic photon and is ex- 
cited to an upper electronic energy level and then proceeds to 
cascade downwards. This is considered populating the energy 
levels from the "top-down", which results in exciting higher H2 
energy levels. In this way, we can use the ratio of certain H2 
lines to discriminate which excitation mechanism is exciting the 
H2 emission. However, the case of UV excitation is sometimes 
degenerate. In the case of high density, the molecules can col- 
lide frequently enough to become thermalized in the lower en- 
ergy levels, which results in line ratios equivalent to the case 



of shock excitation. For the lines to be thermalized, the density 
in the region must be comparable to the critical density of the 
lines, 1.7 x 10^ cm'^ for 1-0 S(l ) and 1.2 x 10^ cm ^ for 2-1 
S(l) (ISternberg & Dalgarndll989l) . 

In order to determine which excitation mechanisms are dom- 
inating which regions, we can use the information from the Bry 
and [Fell] tracers. Specifically, if the H2 is undergoing excita- 
tion by UV photons in a dense environment, this region should 
also be bright in Bry emission if the excitation is by O stars. 
Conversely, if the H2 is being excited by shocks, then the [Fell] 
emission might also be bright in this region. 

In Region 2, coincident with the Bry peak, UV photon ex- 
citation is the dominant mechanism and the density of this re- 
gion is accounting for the lower excitation levels being populated 
leading to the lower H2 line ratio. In Region 6, there is little Bry 
emission. However, the [Fell] emission shows a bright clump. 
Therefore, in this region the emission lines suggest that the H2 
is being excited by shocks. However, in the other regions, there 
is both diffuse [Fell] and diffuse Bry. In this case, we cannot de- 
termine the dominant mechanism by only studying the Bry and 
[Fell] emission. Since H2 can also be excited by slow shocks, 
too slow to destroy dust grains, or by B stars, which are not en- 
ergetic enough to excite Bry, we must use other diagnostics to 
probe the excitation mechanisms. 
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Table 2. Calculated supernova rates for each region and the inte- 
grated galaxy based on [Fell] luminosity. The integrated value is 
integrated over the full map. Supernova rates were calculated 
with the relation betw een [Fell] luminosity and SNrate from 
[Rosenberg eTal] (|2012|) . 



Region 


[Fell] Luminosity 


SNrate 




10" erg S-' 


10"' yr-' 


Region 1 


0.3 


0.04 


Region 2 


0.3 


0.04 


Region 3 


9.8 


1.6 


Region 4 


11.3 


1.8 


Region 5 


17.1 


2.8 


Region 6 


4.1 


0.7 


Region 7 


0.3 


0.04 


Region 8 


0.5 


0.08 


Region 9 


0.4 


0.06 


Integrated 


140.9 


23.5 



One diagnostic is the 2-1 S(l) to 1-0 S(l) transitions at 
2.2478 i-im and 2.1218 jim respectively. This line ratio is com- 
monly used because both lines are relatively bright and have 
similar wavelengths, so they can usually be taken in the same 
spectrum and suffer identical extinctions. Predictions for the ra- 
tio of these two lines due to excitations by shocks range from 0. 1 
to 0.2 dShull & Hollenbach 197 8.), while for photo-excita tion it 
is predicted to be 0.53-0.56 (iBlack & van Dishoecklll987h . 

5.1. Shocks vs Fluorescence 

In Figure|5j we present a map of the 2-1 S(l)/l-0 S(l) ratio with 
the contours of the Bry line map and the rectangular regions 
from Figure[T]overlaid. The 2-1 S(l) line is significantly weaker 
than the 1-0 S(l) transition thus this image has less signal-to- 
noise. We have therefore applied a mask to the map, filtering 
out values lower than 10% of the maximum flux for both line 
transitions. The values of the 2-1 S(l)/l-0 S(l) ratio range con- 
tinuously from to 0.5, showing that we have both thermalized 
and fluorescently excited gas in the galaxy. Both the mean H2 
ratio for each region and the ratio of the integratred 2-1 S(l) and 
1-0 S(l) line fluxes is listed in Table[3] The discrepency in these 
two values is due to the fact that in the left column, we find the 
mean value of the H2 ratio at each pixel for each region, whereas 
in the right column we first integrate the flux in each region for 
the two lines, and then divide the fluxes. The latter method is 
more sensitive to high surface brightness quantities, and shocks 
excited material has higher surface brightness than that of flu- 
orescent excitation (IShuU & Hollenbach| [T978). Therefore, the 
values in the right column are weighted towards shock excita- 
tion, and thus show lower ratios for each region. Averaging the 
ratio on a pixel-by-pixel basis, as done in the left column, does 
not present such a bias. In regions with very high ratios (> 0.3) 
we can confidently conclude that the molecular hydrogen is be- 
ing excited by UV photons from O and B stars. Similarly, there 
are a few regions where the ratio is distinctly lower than the av- 
erage value of ~ 0.2. These regions are region 2, 6, and 7 from 
Figure|5] It is interesting to note that these bright regions all cor- 
relate with an H2 clump or luminosity peak. In these regions we 
cannot determine whether the low H2 line ratio is indicating ex- 
citation by shocks, or a very dense clump being excited by UV 
photons. However, we can place an upper limit on the amount of 
shocked gas. Assuming that all gas with a 2-1 S(l)/l-0 S(l) ratio 
lower than 0.2 is shock excited, we see that 48% of the total 1-0 
S(l) H2 flux may be from "shocked" gas. However, this percent- 



Table 3. Mean H2 2-1 S(l)/l-0 S(l) fine ratios in each region of 
Figure ID (left column) and the ratio of the 2-1 S(l) and 1-0 S(l) 
fluxes (right column). The integrated value is integrated over the 
whole map. 



Region 


/ 2-15(l))\ 
I 1-05(0)) j 


£2-15(1) 
X 1-05(0) 


Region 1 


0.26±0.01 


0.17+0.04 


Region 2 


0.14±0.02 


0.13+0.04 


Region 3 


0.17+0.03 


0.17+0.01 


Region 4 


0.16±0.02 


0.16+0.01 


Region 5 


0.17±0.02 


0.15+0.01 


Region 6 


0.15±0.02 


0.14+0.003 


Region 7 


0.15±0.04 


0.15+0.01 


Region 8 


0.21±0.02 


0.16+0.003 


Region 9 


0.23±0.01 


0.08+0.01 


Integrated 


0.27±0.06 


0.13+0.06 



age is including all regions with a ratio lower than 0.2, including 
regions of high density where the gas is fluorescently excited 
but thermalized. We can at least partially separate this effect by 
filtering out the pixels with the highest Bry flux (> 1 x 10 
erg s"' cm"^), since these regions are dense regions of known 
fluorescent excitation. As seen in the Bry contours in Figure |5] 
this affects Regions 2, 3, and 4. Excluding these regions and any 
pixel with high Bry flux, no more than 29% of all H2 is shock 
excited, and the actual amount may be less. 

5.2. Tracers of Fluorescence 

If strong shocks dominate the excitation of H2 in NGC 253, we 
would expect to see a good correlation between H2 and [Fell], 
which we do not. If fluorescence is the dominant mechanism, 
we would expect the Bry morphology to share more similarities 
with that of the H2. We also see no correlation between H2 and 
Bry emission. However, Bry only traces the most massive and 
young stars, predominately the O stars. In order to excite Bry, 
the star must emitting Lyman continuum photons (A < 912A, 
13.6 eV), yet H2 can be excited by lower energy photons as well 
(A < llOOA, 11.2eV). 

Similarly to H2, polycyclic aromatic hydrocarbons (PAHs) 
have a lower excitation energy than HI (A ~ 42 00A for Nc > 50 
where Nc > 50 is the number of carbon atoms: ^akes & TielensI 
(1994)) and thus are excited not only by O stars but also by 
slightly less massive, B-type stars. Therefore, the PAH emission 
traces the more general fluorescently excited population, includ- 
ing excitation by both O and B stars. By comparing the morphol- 
ogy of the H2 to that of PAH emission, we can determine if it is 
indeed the B stars that are predominately responsible for excit- 
ing the H2 regions, and why neither the Bry, [Fell] or K band 
continuum matches its morphology. 

In Figure |6] we present an ISAAC continuum subtracted 
PAH map at 3.21/im. The left panel is overlaid with H2 1-0 
S(l) contours, the middle panel is overlaid with Bry contours, 
and the right panel is overlaid with [Fell] 1.64 contours. The two 
maps were centred with respect to one another by matching the 
K band continuum peak of the SINFONI maps to the 3.21 con- 
tinuum peak of the ISAAC continuum maps. 

From these images, the morphology of the PAH emission 
matches the H2 emission, better than both the Bry and [Fell]. 
The Bry peak is coincident with secondary PAH peak, mean- 
ing that in this region, there are many O stars producing pho- 
tons E > 13.6ey. However, there is little Bry emission coinci- 
dent with the primary PAH peak, showing that these regions are 
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Fig. 5. The ratio of 2-1 S(l)/l-0 S(l) H2 lines with the Bry line map over plotted in contours. The white rectangles are the same 
regions as defined in Figure[T] The color bar values represent the value of the ratio. The offset from the kinematic center at position 
(a, 5)2000 - (00''47™33M4, -25°17'17".52) is given on each axis. Ratio values where the line flux of either transition is below 10% 
of the maximum flux of that line are masked. 



mainly excited by B stars. The [Fell] peak shows little similar- 
ity to the PAH emission. However, the H2 peaks are coincident 
with both the primary and secondary PAH peaks and shares a 
generally similar morphology in the diffuse PAH emission. It 
is interesting to note that the elongated feature seen in the K 
Continuum, [Fell] and Bry maps in Figure |2] is also present 
in the PAH map, but not in the H2. Since we still have a high 
amount of diffuse emission in this region, it is most likely that 
there is H2 emission originating from this feature, but it is just 
relatively dimmer than the nearby H2 peaks, which therefore 
overshadow it. 

In order to test this coiTelation in a more quantitative sense, 
we can perform a pixel-pixel analysis of the correlation of our 
line maps to that of the PAH 3.21 //m emission. After re -binning 
the ISAAC PAH map to match the pixel size of SlNFONl, we 
plot the flux of H2, Bry and [Fell] as compared to the PAH flux. 
Figure |6] depicts these correlation plots, which show how well 
H2, Bry and [Fell] coiTelate with the PAH flux. The Bry vs PAH 
plot shows an interesting morphology, in which the emission 
seems to be coming from two different areas. First, there is an 
approximately linear correlation for lower Bry fluxes. This cor- 
relation is most likely due to the diffuse Bry emission caused by 
the same photons exciting the PAHs. However, there is a second, 
non-linear, nearly vertical component of extremely high Bry flux 
and medium strength PAH flux. These points coiTespond to the 
location of the super star cluster, which is only the secondary flux 
peak for PAH, but the dominating source of Bry flux. The [Fell] 



also seems to coiTelate with PAH flux with a large, but constant 
spread, suggesting that these tracers originate in the same gen- 
eral region but without a direct physical correlation, emphasized 
by the displacement of the [Fell] peak and the PAH peak. The H2 
also shows a linear correlation with the PAH flux. However, this 
coiTelation has a very large spread at low fluxes and gets contin- 
ually tighter as both H2 and PAH fluxes increase. We can see this 
same trend in the morphologies as well. There is a strong coiTe- 
lation, both spatially and numerically for the peaks of the H2 and 
PAH linemaps, however in the diffuse regions the correlation no 
longer holds. 

The strong, linear correlation between H2 and PAH morphol- 
ogy in the central region of NGC 253 gives added confidence 
that most of the central emission in NGC 253 (Regions 2, 3, and 
4) is fluorescently excited. In addition, the correlation between 
the PAH morphology and H2 morphology indicates that the H2 
in this area, is being excited by the same stars that are exciting 
the PAHs, mainly B-type stars. 

6. PDR and Shock Models 

Although the morphological correlation of H2 to fluorescently 
excited PAHs is a strong indicator that this is the dominant mech- 
anism in NGC 253, we can also compare the H2 line ratios in 
each region to the predicted values from both shock and PDR 
models for a more quantitative analysis. In Figure|5] the ratio is 
lower than predicted for pure UV excitation, although the cor- 
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Fig. 6. Top Panels: PAH 3.28 fim continuum subtracted image from ISAAC. Left panel has H2 1-0 S(l) contours, middle 
panel has Bry contours, right panel has [Fell]i.64 contours. The offset from t he kinematic center a t posit ion (o', 5)2000 = 
(OOMT^SS^.M, -25°17'17".52) is given on each axis. The kinematic center is from iMiiller- Sanchez et al" ] (l2010h . Bottom Panels: 
Pixel-pixel plots of PAH flux vs H2 flux (left), Bry flux (center), and [Fell] flux (right). The pixels of the ISAAC PAH observations 
were rebinned to match those of the SINFONI observations. 



relation to the PAH maps indicates that UV fluorescence is the 
dominant mechanism. In order to study this further, we present a 
diagnostic excitation level diagram. This diagram can be seen in 
Figure [8] and plots the upper energy level (E/k) against the col- 
umn density distribution(A^) divided by degeneracy (g) where: 

gj he Aiv, v', J, J') 

Here, Noi!,s(v, J) is the observed column density of a specific 
upper electronic energy level (v,J), gj is the degeneracy of this 
level, A( v,v',J,J') is the Einstein - A radiative transition probabil- 
ity from IWolniewicz et al.l ([1998), and lots is the observed line 
flux for each transition. The column density is normalized by the 
1-0 S(l) column density divided by its degeneracy (gj). In addi- 
tion, we can plot the best fit Maxwell-Boltzmann distribution to 
the V = 1 transitions. The slope of this line tells us the excitation 
temperature in the region. 

In order to directly compare the observed line ratios in NGC 
253 to model predictions, we can use predicted line intensi- 
ties for shock and PDR models and plot them on this diagnos- 
tic energy level diag ram. Figure [7] presents the results from the 
IShull & HoUenbadhfd 19781) shock models. These models calcu- 
late line intensities for shocks moving at 6, 10, and 14 km/s for 



densities of 10^ 10^ 3-10^ cm . In most of the shock models of 
Figure |7] the v > 1 lines (depicted as triangles) lie on or below 
the thermal distribution indicating subthermal excitation of the 
higher levels. The exceptions to this case are the higher veloc- 
ity (10-14 km/s) shocks in very dense environments (n//=10^). 
However, the v > 1 lines have a shallower slope than the thermal 
distribution as shown by the dotted red line, indicating a higher 
excitation temperature for these transitions. In these high veloc- 
ity, high density shocks, there are two different gas components 
( temperature slopes). 

The bottom 6 panels of Figure 111 displays the predicted line 
intensities from lSternber^ & Dalgarnd (Il989h . who model fluo- 
rescent excitation o f H2 in environments w i th vih - 10^-10^ 
and X - 10^ ^nd iBlack & van Dishoeckl (Il987h who model 
pure fluorescent excitation for a large range of the parameter 
space. We have chosen models 4, 6, 9, 17,19, and 20 spanning 
densities of nn - 10" - lO'* cm"-* and radiation intensities of 
X - 10" - 10**. In contrast to the shock models, the PDR models 
show the V > 1 lines lying high above the thermal distribution 
in almost all cases, since UV fluorescence populates levels from 
the top down, the highest energy levels get populated with higher 
probability compared to the thermalized case. In addition, the 
lines belonging to the v - 2 level share the same excitation tem- 
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Fig. 7. Excitation diagrams of shock and PDR model-predicted H2 line column density distributions as a function of upper energy 
level. The column density is normalized to the 1-0 S(l) H2 transition. Diamonds represent lines from v = 1 while the asterisks 
represent points from v > 1 . The solid line is a best fit of the Boltzmann distribution to t he v = 1 points, the slope represents the 
excitation temperature, Tg^- Predictions of shock models from lShuU & Hollenbachl(ll978h with shock speeds 6 km/s (top row) and 
14 km/s (second row) and densities nn-lO^ cm"^ (left), 10"^ cm"^(middle), and 10^ cm"-' (right). Predictions of PDR models from 
[Black & van Dishoeck(,1987.) with 1„„=10"-10'' and densities nn^lO^ cm ^ 10"* cm ^ and 10^ cm ^ and are labeled with B&vD. 
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perature (same slope) as the v = 1 lines, only offset due to their 
higher upper energy level, unlike the shock models where there 
are two distinct excitati on temperatures. This is the s ame for 
each V level, as shown in lBlack & van Dishoeckl (Il987h . A gen- 
eral feature of the UV-excited models is that the rotational tem- 
peratures (i.e., the excitation temperature determined between 
different rotational levels within the same vibrational level) are 
the same in the various vibrational levels. In contrast, the vibra- 
tional temperature (i.e., the excitation temperature determined 
between different vibrational levels) is higher than the rotational 
temperature. The only exception to this pattern is the top-right 
panel of Fig. 9, showing a model with low radiation field but 
very high density. However, this density is only relevant for very 
dense clumps and the the bulk of the molecular gas in NGC253 
is at much lower densities. 

Now that we have presented the shock and PDR models, we 
can compare them to the observational results seen in Figure [8] 
In these plots, only the (> 2cr) H2 line detections are included. 
Each panel of the figure represents the different regions shown in 
Figure [1] In each region. Figure [8] reveals the thermalized v = 1 
lines of varying excitation temperatures. In addition, in all re- 
gions where v > 1 lines are detected, they lie above the thermal 
distribution. However, the case becomes ambiguous in regions 
where we cannot detect v > 1 lines, such as Regions 8 and 9. It 
is possible that in these regions shocks dominate and therefore 
the higher transitions are not being populated, or that in these 
regions, we lack the signal to noise to detect the lines above 2cr. 

In both Region 1 and Region 3, the kinematic center, we de- 
tect two V = 2 lines and one v = 3 line. Although we only have 
two V = 2 lines, it is clear that they have a similar excitation tem- 
perature (slope) to the thermal distribution and do not represent 
a second hotter gas component. Thus, these lines clearly follow 
a PDR distribution. In addition, the v > 2 lines cannot be excited 
through shocks, so the presence of the 3-2 S(3) line also suggests 
fluorescent excitation. Using these criteria and the shock/FDR 
models as a guide, we can conclude that there is fluorescent ex- 
citation in Regions 1, 3, and 4 and most likely in Regions 2, 
5 and 6. Since the signal-to-noise is too low in Regions 7 and 
8, we cannot determine whether shock excitation or fluorescent 
excitation is the dominant heating mechanism in these regions. 
However, based on the integrated spectrum and Regions 1, 2, 
3, and 4, it is clear that UV excitation in PDRs is the dominant 
mechanism for this galaxy, although shocks may dominate in 
some isolated areas. 

7. Conclusion 

From SINFONI observations, we determine the nature of the ex- 
citation of hot H2 gas in NGC 253 using diagnostic emission 
lines. Specifically, we use the K band continuum as a tracer of 
the older stellar population, Bry as a tracer of high-mass, young 
(O-type) star formation, [Fell] as a tracer of shocks, and PAHs as 
a tracer of slightly lower mass (B-type) star formation. Based on 
the Bry and [Fell] emission line maps, we can interpret the H2 
excitation. We find that in most regions of NGC 253, excitation 
by UV photons in PDRs is the dominant mechanism. There are 
3 regions where shock excitation may dominate, but these are 
small relatively isolated regions. Throughout the entire nuclear 
region of the galaxy, H2 is being fluorescently excited. 

We also present a diagnostic energy level diagram which ro- 
bustly differentiates between shock and fluorescent excitation. 
By comparing the observed column densities of each line to the 
proposed PDR and shock models, it is clear that the NGC 253 
H2 is dominated by fluorescent excitation. A few of the bright H2 



regions may lack the signal-to-noise to resolve the higher energy 
transitions, or it is possible that in these isolated regions, shocks 
are the dominant mechanism. In addition, the PAH emission fol- 
lows the same morphology as the H2, adding further confidence 
that through NGC 253, the gas is undergoing fluoresce. 

We determine that a maximum of 29% of the hot molecu- 
lar gas is excited by shocks. Since we are only sensitive to the 
hottest H2 gas, we are probing the surfaces of molecular clouds. 
Within the molecular clouds, where the bulk of the molecular 
gas is, the gas can have a different excitation mechanism. 
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Fig. 8. Excitation diagrams of observed H2 line column density distribution as a function of upper energy level. The column density 
is normalized to the 1-0 S(l) H2 transition. Diamonds represent lines from v = 1 while the asterisks represent points from v > 1. 
The solid line is a best fit of the Boltzmann distribution to the v = 1 points, the slope represents the excitation temperature, Tg^- 
Each panel represents a different region, corresponding to the regions in Figure|5] The error bars represent a standard deviation from 
the column density. 
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